Spin current and magneto-electric effect in non-collinear magnets 
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A new microscopic mechanism of the magneto-electric (ME) effect based on the spin supercurrent 
is theoretically presented for non-collinear magnets. The close analogy between the superconductors 
(charge current) and magnets (spin current) is drawn to derive the distribution of the spin super- 
current and the resultant electric polarization. Application to the spiral spin structure is discussed. 
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The interplay between the magnetism and ferroelec- 
tricity is an old issue since the firstprediction of the 
magnetoelectric (ME) effect by Curie Later the phe- 
nomenological theory of ME effect was developed by 
Landau 2] and Dzyaloshinskii^. There the symmetry 
consideration is essential to classify the all possible ME 
tensors depending on the magnetic point group. Es- 
pecially the time-reversal (T) and spatial inversion (1) 
are the key symmetries to the ME effect. For example, 
the linear ME effect corresponding to the term aijMiPj 
(P:polarization, A/ magnetization) is allowed only when 
both T— and /—symmetries are broken. However, on the 
other hand, the microscopic quantum theory of ME effect 
has not yet been fully developed although several sce- 
narios for particular materials such as Cr2 03 have been 
proposed 4, 5]. For this particular material, the change of 
the anisotropy energy, exchange, and y-value due to the 
electric field have been proposed for the origin of the par- 
allel ME effect. The transverse ME effect, on the other 
hand, can not be explained by these mechanisms. For a 
spiral spin magnet such as ZnCr2Se4, the electric field in- 
duced Dzyaloshinskii-Moriya (DM) interaction 0,0 has 
been proposed for the transverse ME effect 0. Another 
well-studied ME material is (Ga,Fe)03, where the ME co- 
efficient is an order of magnitude larger than in Cr203, 
which appears mostly in the anti-symmetric transverse 
channel. This ME effect was attributed to the presence 
of toroidal moment T = ^ Eq x .Sq [9]. 

We will argue that ME effect and spin current are di- 
rectly related. The spin current js has attracted revived 
interests recently in the context of spintronics in semi- 
conductors. In contrast to the charge current, it is T- 
even since the spin polarization is also reversed together 
with the direction. Therefore from the viewpoint of pure 
symmetry, js belongs to the same class as the electric 
polarization P, and it is natural to expect the coupling 
between these two. In fact, the electric field induced 
dissipationless spin current driven by the Berry phase 
curvature has been proposed for the semiconductors |lfll | . 
In magnets, the spin current is analogous to the super- 
fluid current, i.e., spin supercurrent, associated with the 



spin rigidity In this paper, we present a new mech- 
anism for the ME effect due to the spin supercurrent 
with the non-collinear spin structure such as the spiral 
state^2 The idea is that the spin supercurrent is 
induced between the two spins with generic non-parallel 
configurations, which induces the electric polarization. 
This is the dual effect to DM interaction[a, and/or 
Aharonov- Casher effect Q. The vector potential cou- 
pled to the spin current is the vector product of electric 
field/electric polarization and the direction of the bond 
connecting the two spins [l^ . This allows one to develop a 
microscopic theory of the electric field induced DM inter- 
action (^]. Applications to the spiral magnetic structure 
are also discussed. 

Let us first review the analogy of the magnetically or- 
dered state to the superconducting state. The key to 
understand the superconductivity is the canonical con- 
jugate relation between the charge n and the Josephson 
phase ip, i.e., [ni,tpj] = iSij, where i,j are the indices 
for the site. Similar relation exists in the quantum spin 
operators as [S'f , 6j] — iSij, where 5*^ is the z-component 
of the spin operator while 9 is the angle of the vector 
(5^, S^) measured from the cc-axis. This makes the link 
between XY spin model and that of superconductivity. 
The Hamiltonian for XY model reads 

HxY^-J2^iS+S- + Srs;) (1) 

<ij> 

and the spin supercurrent jfj defined so as to satisfy 
dSf/dt = il/th)[S^,H] = - J|. is given by 

j^^=tJ^,,iS+Sr -SrS+) (2) 

Putting iS^,S^) = S{cosej,sin0j), We obtain 

jt,=J±^,S^smi^,-^,) (3) 

where J±S^ corresponds to the spin stiffness, i.e., rigidity. 
Equation (3) is analogous to the Josephson equation. 

To go further with this analogy, the next question 
is "what is the vector potential Ag coupled to the 
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spin supercurrent ?". The answer to this question can 
be found in the Aharonov-Casher (AC) effect and 
Dzyaloshinskii-Moriya (DM) interaction [^Q- The con- 
ventional DM interaction[^ is given by 



H 



DM 



<»i> 



A, • [S^ X S,). 



(4) 



When the DM vector Dij ~ DijCz, the total Hamiltonian 
^^totai = HxY + Hdm with HxY in eq. (1) is written as 



total — - ^ — ^ (e 
<i]> 



S+SJ +^^^^StS-) (5) 



where J^ij e'"**^ = J^_ij + iDij . Therefore the DM vec- 
tor D acts as the vector potential or gauge field to the 
spin current. It is well known that the DM interaction 
exists only when the inversion symmetry is broken at the 
middle point between the two spins. Therefore when the 
crystal structure has the inversion symmetry, the exter- 
nal electric field E induces the DM interaction. Namely 
Di.j cx E X e^j, where is the unit vector connecting the 
two sites i and j. This form is identical to the Aharanov- 
Casher (AC) effect, where the Lorentz transformation of 
the electric field induces the magnetic field in the moving 
frame which interacts with the spin moment. However 
the magnitude of the coupling constant for AC effect is 
extremely small in vacuum since it contains the rest mass 
of the electron mc^ = 5 x IQ^eV in the denominator. The 
situation is different for the DM interaction in solids, i.e., 
the electrons are trapped in the strong potential of the 
atoms with large momentum distribution leading to the 
enhanced spin-orbit interaction. Therefore the gauge po- 
tential Aij could be (a fraction) of the order of unity, e.g. 



A,. 



2tt as seen below. 



To illustrate this, consider the electron energy levels in 
the ligand field of 3d-transition metal 0|. In the octa- 
hedral ligand field, the d-orbitals are split into Cg orbitals 
and t2g orbitals. The t2g orbitals, i.e., dxy^ dyz, and dzx, 
have energy lower than Cg orbitals. If we take account 
of the spin degree of freedom, there is 6-fold degeneracy 
in t2g energy level. Due to the on-site spin-orbit interac- 
tion, however, this degeneracy is lifted and we have two 
groups of spin-orbit coupled states, labeled Tr and Fg. 
The 2-fold degenerate states, i.e., Ty, are given by 



the above two states alone. However, our method is valid 
for more general cases and one can easily generalize it to 
any other spin-orbit strongly coupled situation. 

We consider the case where the inversion symmetry 
exists at the middle point of the two magnetic ions, and 
the generic non-collinear magnetic ordering is realized 
by the competing exchange interactions J's and/or by 
the symmetry breaking due to the spin-orbit interaction. 
Here the magnetic moment at j-th site points to the unit 
vector ej = (cos sin 6j , sin sin 9j , cos 9j ) . The mean 
field Hamiltonian applied to the Hubbard model is given 
by ( we take the unit where h ^ 1 hereafter): H — 
~UJ2j ■ ^jj where U is energy of Coulomb repulsion. 
For each site j, we restrict the Hilbert space to the 2- 
dimensional one spanned by the above two states, and 
the effective Hamiltonian is reduced to the 2x2 matrix 



U 
3" 



smfe 
cos^ 



(8) 



We diagonalize this Hamiltonian matrix to obtain eigen- 
states |P), \AP) as 



IP) 



\AP) 



sm - a) 
2' ' 



cos -I a) 



e^-^ cos ■ 



e^^sin-|6). 



(9) 



Here \P) and \AP) means the spin state parallel and 
anti-parallel to the unit vector e, and the corresponding 
eigenvalues are — -j and respectively. For conve- 

nience, we define the coefficients A^"' and B^"' and ab- 
breviate the above two states as, \P) — X^itr ^*'^Mio'), 
\AP) = Y^iaB^'^ldia), where i = xy,yz,zx,(j =],l. 

From now on, we focus on the three atom model as 
shown in Fig.l, which represents the bond between the 
two transition metal ions Ml, M2 through the oxygen 
atom O. We take the hole picture below, where the oxy- 
gen orbitals are empty. We assume the generic case of ei 
and 6*2 including the non-collinear configuration. Each 
site has two states, i.e., \P) and \AP), mentioned above. 
So we define \P)j and \AP)j {j = 1, 2) corresponding to 
the magnetic order on each site. Because of the existence 
of the oxygen atom, there are hopping processes between 
the M site and the O site. The transfer integrals between 
the d- and p-orbitals can be found in the Slater-Koster 
tables [tR Il8j|. and the hopping Hamiltonian is given as 
follows: 



|a) = -^i\dxy^^) + \dyz,i) + i\dzx,i)), (6) 



tit + tit 



and 



\b) = ^{\d.y,i) - \dyzA) + «|4.,t)), (7) 



respectively, where the quantization axis of spin is taken 
to be the z axis. For the sake of simplicity, we consider 
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d-orbitals 



FIG. 1: The cluster model with two transition metal ions 
Ml, M2 with the oxygen atom O between them. With the 
non-coUinear spin directions ei and 62, there arises the spin 
current js cx e 1 x 62 between Ml and M2. Here the direction 
of the vector js (denoted by the green arrow) is that of the 
spin polarization carried by the spin current. The direction 
of the electric polarization P is given by P cx ei2 x js where 
ei2 is the unit vector connecting Ml and M2. 



where V(> 0) is the transfer integral and the superscript 
j denotes the corresponding site number. Let us now 
treat the above Hamiltonian perturbatively in V. The 
eight bases we must prepare for this problem are \P)j, 
\^P)j, (i = 1,2), andpi,^, (i = y, z, cr =t, i). Using the 
second-order perturbation theory, the four lowest lying 
states and corresponding perturbed energies are obtained 
as follows: 



ID 



ae 



V2\c 



with El = -|^(l + |a| 



|2> = - 



with E2 = ^^^{^ ^ |«|), and two other higher energy 
states. Here A(> 0) is the energy difference between 
the p-orbitals and \P)j, Acj) = (pi — 4>2, and we have in- 
troduced the complex number a = cos ^ cos ^e~'^^ + 

sin sin ^e+'~ . Before calculating the expected value 
of the polarization, it is useful to note that only the fol- 
lowing matrix elements are non-zero from the shapes of 



d- and p-orbitals; 

/= / d^fdl^){r)yp,{r), (i = l,2), 



and its cyclic permutations. The integral / is approxi- 
mately estimated as / = ^Zq^^^^ Zu'^^^i^ + ^)"^ao , 
where ag is Bohr radius and Zq/Zm is the atomic num- 
ber of 0/M. We can easily check the above results by 
expanding wavefunctions in terms of lattice constant a. 
So let us now calculate the expected value of polarization 
in the following two cases. 

Double- exchange intera,ction\l^ l2(Tj|. — First, we con- 
sider the situation where only one hole is present. In this 
case, this hole is put into the ground state, determined 
by the above second order perturbation theory, and the 
expected value of polarization, (l|er|l)/(l|l), is given by 



P 



eV j. ei2 X (ei x 62) 
'3A Icos^l 



(10) 



where 612 is the unit vector parallel to the direction of 
the bond from site Ml to site M2, and 9x2 is the an- 
gle between the two vectors ei and 6*2, i.e., e 1 • 6*2 = 
cos 6112 . Spin current Jg is approximately given by Jg ~ 
(y^/A)ao(ei X 62)/ cos(6'i2/2), and eq.(lO) can be rewrit- 
ten as P ~ (e/y)ei2 x js. Therefore the spin current is 
essential to the electric polarization. 

Superexchange interaction 2\\. — Next we consider the 
case of two holes. From a viewpoint of Hartree-Fock ap- 
proximation, two holes are put into the ground state |1) 
and the second low- lying state |2) , and the expected value 
of the polarization is given by the following form similar 
to (10): 

^=-f (^m2x(eixe2). (11) 

In this case, we must pay attention to the difference of the 
normalization factor between the two perturbed states, 
i.e, (l|l>-i - 1 - 1(^)2(1 - \a\) and (2|2)-^ = 1 - 
|(^)^(1 + |q;|). Of course there are many other terms but 
this term mainly contributes. As the above result is order 
of (F/A)*^, one may think that the polarization is too 
small to observe. By using the superexchange interaction 
J, i.e., J = y^/(C/A^), however, (11) can be rewritten as 

4e J U 

P = -y y A^e^i2 X (ei X 62)- (12) 

Again this equation can be interpreted in temrs of the 
spin current js as P ^ {eU/VA)ei2 x js- Therefore the 
magnitude is not too small on the ground that the ratio 
of [/ to A is by no means small practically. 

For more general models, the magnitude of the elec- 
tric polarization induced by the spin current would differ 
from those obtained in the present model, and depends 
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(b) 




0^ 



V7 If ~ 

J, 



(d) 



FIG. 2; (a) Generic spiral spin configuration, (b-d) Some 
of the specific configurations wliere tfie geometrical relation 
among spins(black arrows), spin current(green arrows), and 
electric polarization (blue arrows) are shown. 



on the details of the electronic level structure. Most likely 
it will be smaller than in the present case, especially for 
the eg systems where the spin-orbit interaction is smaller. 
However the geometrical relation between the spin cur- 
rent and electric polarization remains unchanged. 

Applications to the spiral magnets — Now we turn to 
the discussion on realistic non-coUinear magnets. One 
of the typical examples is the spiral structure where 
the direction of the spin rotates along the wavevec- 
tor q. Figure 2(a) shows the most generic spiral spin 
configuration where the spiral axis is along the x-axis, 
the cone axis direction have the angle a measured 
from the z-axis, and the cone angle is /3. The an- 
gle 7j for the j-th spin direction is measured from 
the zx-plane, and the undistorted spiral means 7^ = 
QJ + 7o- The spin at site j can be written as Sj = 
^(cos P sin a + sin f3 cos a cos 7j , sin /? sin , cos f3 cos a — 
sin /3 sin a cos 7j ) . With this configuration one can 
easily calculate the spin current jsj+1/2 and resul- 
tant electric polarization Pj+1/2 at each link connect- 
ing j and j + 1 as {Pj+i/2)x = 0, {Pj+i/2)y ^ 
— cos (3 sin/3 sina[sin7j+i— sin7j]— sin^ /3 cos a sin(7j+i — 
7j), {Pj+i/2)z ~ cos/3sin/3[cos7j+i - cos7-,]. There- 
fore only the y-component of the uniform electric po- 
larization P = '^jPj+1/2 is nonzero and is given by 
Py oc J2j sin^ P C0SQ!sin(7j4.i - 7^). 

Figure 2 (b-d) shows some typical cases and their spin 
current and electric polarization. When the spins are 
within the xy-plane as in Fig. 2(b), i.e., a = 0, /3 = 7r/2, 



the spin current js is along the z-direction, and the 
electric polarization P is along the y-direction for each 
site. Therefore the total uniform polarization is finite 
along the y-direction. Note that even when the spiral 
wavenumber q is incommensurate with the lattice peri- 
odicity, the uniform polarization, i.e., the ferroelectric- 
ity, is realized. When the spins are in the yz-plane (Fig. 
2(c)) , i.e., a = 7r/2, P = tt/2, both ejj+i and js are in 
the x-direction and their vector product is zero. There- 
fore we do not expect any electric polarization induced 
in this case. Figure 2(d) shows the "conical" spin struc- 
ture, where the finite x-component of the spin is induced 
starting from the structure in Fig. 2(c), i.e., a = 7r/2, 
< P < ■n 12. In this case, the finite S^, component at 
each site produces the rotating polarization, but these 
cancel out to zero uniform electric polarization. 

Now consider the effect of the external electric field 
Ey along the y-direction on the generic spiral configura- 
tion. This field induces uniform magnetization (per site), 
estimated as fh^^ oc i?y sin^ /Jcos/Jsinacosasing, and 
TO2 (X sin^ /3cos/3sin(7(a + 6cos^ a) (a, 6: constants). 
This dependence is consistent with the group theoretical 
consideration and the experiments on ZnCr2Sc4 ( Figures 
6 and 7 of ref. 8| ). 
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